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We analyze the problem of quantifying entanglement in pure and mixed states of fermionic sys¬ 
tems with fixed number parity yet not necessarily fixed particle number. The “mode entanglement” 
between one single-particle level and its orthogonal complement is hrst considered, and an entan¬ 
glement entropy for such a partition of a particular basis of the single-particle Hilbert space H is 
defined. The sum over all single-particle modes of this entropy is introduced as a measure of the 
total entanglement of the system with respect to the chosen basis and it is shown that its minimum 
over all bases of "H is a function of the one-body density matrix. Furthermore, we show that if 
minimization is extended to all bases related through a Bogoliubov transformation, then the entan¬ 
glement entropy is a function of the generalized one-body density matrix. These results are then 
used to quantify entanglement in fermion systems with four single-particle levels. For general pure 
states of such a system a closed expression for the fermionic concurrence is derived, which generalizes 
the Slater correlation measure dehned in [J. Schliemann et al, Phys. Rev. A 64, 022303 (2001)], 
implying that “particle entanglement” may be seen as minimum “mode entanglement”. It is also 
shown that the entanglement entropy dehned before is related to this concurrence by an expression 
analogous to that of the two-qubit case. For mixed states of this system the convex roof extension of 
the previous concurrence and entanglement entropy are evaluated analytically, extending the results 
of previous ref. to general states. 
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I. INTRODUCTION 

Quantum entanglement is not only one of the key fea¬ 
tures of quantum mechanics but also an essential resource 
in quantum information processing [l| . It plays a central 
role in quantum teleportation Q and quantum computa¬ 
tion 0 . Consequently, the understanding and quantifica¬ 
tion of this resource has become a fundamental problem 
in quantum information theory [^. It has also provided 
deep insights into the structure of correlations and quan¬ 
tum phase transitions in many-body systems im. 

If I'Fab) S FLa ^FLb is a pure state of a composite 
quantum system, its entanglement is quantified by the 
entanglement entropy S{pa) = S'(TrB|4'AB)(4'TB|) = 
S{pb), where S{p) = —Trplog 2 P is the von Neumann 
entropy. It is then seen that the notion of entanglement 
in such systems relies on the tensor product structure of 
its state-space Q. In fermionic systems, however, the 
situation is less clear since the state-space has no longer 
this structure due to indistinguishability. 

When generalizing the notion of entanglement to sys¬ 
tems of indistinguishable particles imi mainly two dif¬ 
ferent approaches have been taken: Mode entanglement 
fl3 . a nd qua ntum correlations/particle en¬ 
tanglement [OMlll. [isHTM [2ll - l^ . In the first case the 
parties share different modes of a given basis of the 
single-particle Hilbert space. Therefore, mode entangle¬ 
ment of a system does not remain invariant with re¬ 
spect to unitary transformations in the single-particle 
(sp) space. The second approach looks for correlations 
between particles and beyond antisymmetrization. In 
i[iil a fermionic analog of the Schmidt decomposition 
and Schmidt number was introduced to quantify entan¬ 
glement in two-fermion systems, and also a fermionic 


“concurrence” was defined. While these measures of en¬ 
tanglement remain invariant under unitary transforma¬ 
tions in the sp space, they are restricted to states with 
a fixed particle number, which is not the general case 
in fermionic systems. The same problem arises in [l5l| . 
where in order to share particles between parties it is nec¬ 
essary to project the original state onto subspaces with 
definite particle number. 

In this paper we first consider pure states of fermionic 
systems within a grand-canonical context, so the particle 
number is not necessarily fixed. Fermionic states with no 
fixed number of fermions arise, for instance, when con¬ 
sidering the vacuum of quasiparticles defined through a 
Bogoliubov transformation [^, , as well as by sim¬ 

ply applying particle-hole transformations, such that the 
state is viewed as a vacuum of certain fermion operators 
plus particle-hole excitations. The fermion number par¬ 
ity of these states is nonetheless fixed, in agreement with 
fermionic super-selection rules psj j. The entanglement 
between a single fermionic mode and the remaining sp 
orthogonal space in such states is first considered, and 
an entanglement entropy is defined in order to quantify 
these correlations. We then propose the sum over single¬ 
particle modes of this entropy as a measure of the total 
mode entanglement associated with the chosen sp ba¬ 
sis, and show that its minimum over all single particle 
bases depends only on the eigenvalues of the one-body 
density matrix = (c|ci), being therefore invariant un¬ 
der sp transformations. Furthermore, it is shown that if 
the minimization is extended to all quasiparticle bases, 
i.e., bases related through Bogoliubov transformations, 
the minimum entanglement entropy is just the von Neu¬ 
mann entropy of the generalized one-body density ma¬ 
trix which contains in addition the pair creation 
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and annihilation contractions (cjc|^) and (cjCi). Its con¬ 
vex roof extension for mixed states is also introduced. 
This quantity allows to rigorously identify mixed states 
which cannot be written as convex mixtures of Slater de¬ 
terminants or quasiparticle vacua (or in general fermionic 
gaussian states [^,| 23 )) like thermal states of interacting 
fermion systems at sufficiently low temperatures, quan¬ 
tifying their quantum correlations. 

We then focus on fermionic systems with sp space di¬ 
mension 4. For general states it is shown that the mini¬ 
mum over all quasiparticle bases of the entanglement en¬ 
tropy can be written in terms of a fermionic analog of the 
concurrence [1^ [H, [2^ , that reduces to the Slater corre¬ 
lation measure defined in iini for two-fermion states. 
Its convex roof extension for mixed states is also eval¬ 
uated analytically, extending explicitly the results of [p 
to arbitrary mixed states with fixed number parity [22j |. 
This allows to evaluate in closed form the convex-roof ex¬ 
tension of the previous entanglement entropy. A simple 
illustrative example is provided. 

II. FORMALISM 

A. Single level entanglement entropy 

We start by considering a pure state |4') of a fermion 
system with an n-dimensional single-particle (sp) Hilbert 
space H. The system is described by a set of fermion 
annihilation and creation operators {cj,cj,} satisfying 

CiCj CjCi = 0 , CiCj -|- CjCi = Sij , ( 1 ) 

such that {c||0), j = 1,... ,n} is an orthonormal basis 
of sp states (|0) denotes the vacuum of the operators Cj). 
We will work within a general grand-canonical context, 
in which |'I') is not necessarily a state with a definite 
value of the fermion number N = be, 

for instance, a vacuum of qusiparticle operators a,^, re¬ 
lated with the Cj ’s through a Bogoliubov transformation 
[^ . In this case, it is a sum of pure states with different 
fermion numbers (see Appendix), yet having all the same 
number parity 

P = exp[f7r^c|cj], (2) 
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such that P\'^) = il^*). Let us also recall that the ele¬ 
mentary particle-hole Bogoliubov transformation 

Cj —>■ cj, cj —>■ Cj , (3) 

leaves the anticommutation relations unchanged, so that 
formally, it is a matter of choice whether one considers 
the particles or the holes as the “true” fermions. We will 
take this basic symmetry into account in all the following 
correlation measures, such that they all remain invariant 
under the previous transformation. We will just assume 


that all pure states involved have a definite number par¬ 
ity [1^, which implies (cj) = 0 and also (O) = 0 for 
any operator O which is a product of an odd number of 
fermion operators Cj , cj. 

We now consider a partition {A,B) of TL, where A 
denotes the single mode or “level” j and B the remaining 
orthogonal sp space. Eq. © implies that the operators 

Uj = cjcj, Uj = Cjc] , Bj -k Ilj = 1, (4) 

constitute a basic set of orthogonal projectors, defining 
a standard projective measurement on the level j. Ac¬ 
cordingly, we may decompose any state I'k) as 

1^) =cjc,|4')+c,cj|4') (5) 

= >) ( 6 ) 

where the first (second) term in ([5]) selects the compo¬ 
nent of |4') where the state j is occupied (empty) and 
I'^j) = 1%) = ^Cjcjl^-) are the corre¬ 

sponding normalized states. Here pj (pj) is the probabil¬ 
ity of finding the level j occupied (empty) in |'L): 

P3 = (^Icjffi |4'), p-j = (4'|cjcj|4') = 1-Pj. (7) 

For an operator Oa depending just on Cj, cj and Ob de¬ 
pending just on the complementary set {cfc,cj,,fc j}, 
we then obtain, assuming PIT) = ±|T), 

(T|0^(s) |T) = p, (T, |0^(s) IT,) + |Tj) 

= Ba(b)Pa{b)Oa{b) , ( 8 ) 

where pA = p,cj|0)(0|c,- + pj|0)(0| and pB = 

р, c,jT, )(T, |cj-|-pj|Tj)(Tj| represent reduced density op¬ 
erators for systems A and B respectively. 

The entanglement between A and B can then be 
quantihed by the entropy of the elementary distribution 

{pbP 3 = ^-Pj}' 

S{pa) = S{pb) = -pj\og2Pj - (1 -p,)log2(l -p,)(9) 
= HPj) , (10) 

where S{p) = —Trplog2/5 is the von Neumann entropy 
and/i(p) = —plog 2 P—(1—p) log 2 (p) (0 < h{p) < 1). This 
entropy remains obviously invariant after a particle-hole 
transformation ([21). For a pure state |T), Eq. ([5|) vanishes 
if and only if (iff) |T) is separable with respect to this 
level, i.e., iff the level j is either occupied {pj = 1) or 
empty {pj = 0) in |T), such that |T) = cjc,jT) or |T) = 

с, cj|T) respectively. Its maximum value 1 is attained for 

p, = 1/2. 

B. One-body entanglement entropy 

The sum 

So = J2Hpj), (11) 

j 
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is a measure of the entanglement associated with the sp 
basis determined by the operators cj. Eq. m vanishes 
iff each level j of this basis is disentangled from its com¬ 
plementary sp space, i.e., iff each level is either occupied 
or empty in such that is a Slater determinant in 
this basis: |'I') = .. -cJ^lO) for some subset of levels 

{jl; ■ ■ • ; Jm}- 

Eq. m depends on the choice of sp basis, i.e., on the 
choice of fermion operators c = (ci ... ,c„)^. We now 
consider the minimum of (ITTI) over all sp bases of TL, i.e., 

S^P =MinS'e^, (12) 

c' 

where Sc' = P'j = (^|c'Jc'and c' = 

(cj,..., an arbitrary set of fermion operators related 
with the Cj’s through a unitary transformation: 

c' = U^c, (13) 

with U a n X n unitary matrix (such that the fermionic 
relations 0 are preserved). Eq. (fT^ vanishes iff \'i') is 
a Slater determinant i.e., j^*) = c'\^ .. .c'j,^|0) for some 
operators c). of the form (fT^ . Hence, S^'^ = 0 iff there 
is a sp basis where every level is disentangled from its 
complementary sp space. 

Defining the sp density matrix p®P = 1 — (ccf) (with 
(O) = (5'|0|4')), of elements 

Pfj = {c]ci ), (14) 

it is seen that the minimum m is reached for those 
operators c' which diagonalize p®P, i.e., satisfying 

(c'Ic() = (C/tp^P[/),fc=p;4i, (15) 

with p'j, the eigenvalues of p’^P. 

Proof: Eqs. (II31)-(ll5l) imply that pj = pfj = X);, Pfklp'k- 
Hence, concavity of the function h{p) entails h{pj) > 

J2j.k\Ujk\HPk) = J2k Hp'k)^ with equality reached iff the 
Pj's are already the eigenvalues of p®P. □ 

The minimum value (1121) can then be expressed as 

= = (16) 
k 


being now apparent that ^^p vanishes iff the eigenvalues 
p'j, are either 0 or 1, i.e., iff (p®P)^ = p‘^p, a condition 
ensuring that |'I') is a Slater determinant 241. 

Eq. (1161) has in addition the obvious meaning of quan¬ 
tifying how mixed or “hot” is |4') with respect to the set 
of all one-body operators of the form 


C. Generalized one-body entanglement entropy 

A quasiparticle vacuum, like for instance a superfluid 
or superconducting state in the BCS approximation [^ . 
will lead to S'®p > 0, since p®P will be mixed, i.e., it will 
have eigenvalues distinct from 0 or 1 (see Appendix). If 
fermion quasiparticles are to be allowed, we can extend 
the minimization in (1121) to all single quasiparticle basis, 
i.e.. 


S^i'P = MinS'a, (18) 

a 

where Sa = h{{aXaf)) and a denotes a set of fermion 
operators linearly related to the original operators Cj , 
cj through a general Bogoliubov transformation [^ : 

Oi, = 'y ( UjiiCj -\- ■ ( 19 ) 
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Eq. (fTOl) can be written as 

where the 2n x 2n matrix W should be unitary (i.e. 
Dt/t -b VV^ = 1, + VU'^ = 0) in order that the 

operators Ojy, fulfill the fermionic anticommutation re¬ 
lations (1). 

One should then consider the extended 2n x 2n density 
matrix 

p*"I-((,■;)(-• ■=))=(4 i_v).(2i) 

where k is an n x n antisymmetric matrix containing the 
pair annihilation averages 

tvij — (,CjCiJ , (^^) 

with -Kij = (c]cj) and (1 - = (cjcj). Eq. (I2T|) is 

a hermitic matrix which can always be diagonalized by a 
suitable transformation (I20L such that 

1-((C,)(» a')) = WV-W = ( ( i4), 

(23) 

with = fi^Sj^i, and /j,, 1 — fv the eigenvalues of p'^^P 
(which always come in pairs (/,y, 1 — /i,), with fi, S [0,1]), 
entailing 


0 = Ylo^jc\cj, (17) 

i,3 

since their averages are completely determined just by 
psP: (TjOI'l') = tr p^Pq. Accordingly, S'®p remains in¬ 
variant under one-body unitary transformations |'I') —>• 
exp(—iO)|'I'), with O any hermitian one-body operator 
of form (HZl, since they lead to a unitary transformation 
of p®P (p®P —>■ C/p®P17t, with U = e~'°) and hence do not 
affect its eigenvalues. 


(afajf) — djii/fi,, {ajiaf) — 0 . (^ 4 ) 

It can then be easily shown that the minimum (1181) is 
^-qsp = _ ^ log^ fu + {l- U) log2(l - U) (25) 

U 

= -tr' p'l'P log 2 p'l^P . (26) 

where tr' denotes the trace in the extended sp space. 
Proof: Since bothpj = (cjc^) and 1—pj are the diagonal 
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elements of denoting with qj and the full set 

of diagonal elements and eigenvalues of we obtain 
qj = hence, due to the concavity of 

f{p) = -plogaP, S'c = > J2j,^\^jJfM = 

= □. 

Eq. (1^51) vanishes iff fi, is either 0 or 1 for all i', i.e., iff 
is a particle or quasiparticle Slater determinant. By 
an elementary particle-hole transformation we can always 
change such state to a quasiparticle vacuum, so that we 
can say = 0 iff 14') is a quasiparticle vacuum. In 
other words, 5"^®^ = 0 iff there is a single quasiparticle 
basis where every level is disentangled from the rest. 

Eq. (fT51) also measures the mixedness of I'h) with re¬ 
spect to the set of all generalized one-body operators, of 
the form 


= X] + W1°CiCj + o°]c\c]) - itr o“ (27) 




3^ 

^11 

o20 


^02 


m 


i.e., general quadratic functions of c, <X (the constant 
term in (|27)) is just added for convenience), since their 
averages are completely determined by 


(^|0|^) = tr[p"P o“-io“ + i(Ko20-Ko°2)] = itr'p'i'PO. 

(29) 

The present scheme allows then to properly treat states 
which do not have a definite fermion number and lead 
to non-zero contractions (ciCj). The whole formalism 
becomes then strictly invariant under arbitrary particle 
hole transformations (jS]) applied to some subset of levels, 
which will move elements from p®P to k and viceversa, 
but which will not alter the spectrum of p'l^P. The latter 
remains actually invariant under arbitrary quasiparticle 
unitary transformations |'I') —>• exp[—zO]jd^), where O is 
an hermitian generalized one-body operator of the form 
dSll), since they just lead to a unitary transformation of 
p'l'P, i.e., p'l'P Wp‘^"P>V^ with W = e"*®. 

We notice that a transformation a^, O al obviously 
changes /i/ O 1 — /^, so that there is no unique way to se¬ 
lect which of the eigenvalues of p'^®P will be the /i/’s or the 
1 — /'s. One can choose the /j/’s as the lowest eigenval¬ 
ues (such that 14') becomes a quasiparticle vacuum when 
^qsp _ jg g^jgp possible to set Det[I/] ^ 0, which 

ensures that the vacuum of the has the same number 
parity as |0) (Eq. (lAlll l. These choices do not affect the 
entropy 5 ''^®p. We also remark that the maximally en¬ 
tangled state, i.e., that with maximum 5 "^®p, corresponds 
to the exceptional case f^, = 1/2 V v, where = n 
and p‘i®P = hn/^ becomes proportional to the identity 
matrix, remaining then invariant under any Bogoliubov 
transformation. 


D. Generalized entropic inequalities and quadratic 
entropy 

From their definitions, it follows that the entropies 
(HH), (HU) and (1^ satisfy the inequality chain 

. (30) 

Eq. (1301) actually holds for more general entropic forms. 

If p®P = ^ ^ ^ is the extended p®P and 

the diagonal of p®P, we obtain, with the same previous 
arguments, 

Sf{p7)>Sf{p^P)>Sf{p^^P), (31) 

where 

Sf{p) = tTf{p), (32) 

with / : [0,1] —>■ K. a strictly concave function satisfying 

/(O ) = /(I) = 0, represents a generalized entropic form 
Moreover, these matrices fulfill the majorization 
relation [2^ 

p7 -< -< . (33) 

where p ^ p' means here Y7i=i — EUK for J = 
1,..., 2n — 1, with Xi, A' the eigenvalues of p and p' sorted 
in decreasing order, since the sorted set of diagonal ele¬ 
ments in an orthonormal basis of an hermitian operator 
are always majorized by the sorted set of its eigenvalues 
. Eq. (15^ allows to extend m to any Schur-concave 
function [29| of the extended density matrices. 

A particularly useful example, which will play an im¬ 
portant role in the next section, is the quadratic en¬ 
tropy S 2 ip) (also denoted as linear entropy), obtained 
for f{p) = 2p{l-p): 

52(p‘ 1"P) = 2tr'[p'i"P(l-p‘i"P)] 

= 4tr[p®P(l — p®P) — Kf^t] (34) 

= 4^/,(l-/,), (35) 

ly> 

where the factor 2 has been chosen such that its maxi¬ 
mum value for a single level is 1. Unlike the von Neum- 
man entropy (1261) . S 2 {p^^P) can be evaluated just by 
taking the trace in (I34L without explicit knowledge of 
the eigenvalues /y of p^^P. Yet, like it is non¬ 

negative and vanishes iff |'I') is a quasiparticle vacuum 
or Slater Determinant. Eq. (EH implies in particular 
EjPji^-Pj) > EkP'ki^-p'k) > - !v)- 

E. Mixed states 

Let us now consider mixed fermion states, assumed 
as convex mixtures of pure states with definite number 
parity, i.e., 

P = '^qi\'^i){'^i\, (36) 
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where qt > 0, '^^qi = 1 and -P|'f'i) = il'I'i), such that 
\p, P] = 0. We can define an entanglement measure for 
these mixed states in a way similar to the entanglement 
of formation Bmi, through the convex roof extension 
of 


^ Mm ^ ^ , (37) 


where p = q'i\^'i){^'i\, q'i > 0, and the minimization is 

over all decompositions of p as convex mixtures of pure 
states, assumed again of definite number parity. Eq. dSil) 
vanishes iff p is a convex mixture of particle or quasipar¬ 
ticle Slater determinants, i.e., of suitable quasiparticle 
vacua, and reduces to for pure states. This quantity 
will be evaluated exactly in the particular system of the 
next section. 

As a general application of let us consider an 

interacting fermion system at finite temperature T. For 
attractive two-body couplings, the static path approxi¬ 
mation [H, will lead to a classically correlated den¬ 
sity operator psPA, which is a convex mixture of (non¬ 
commuting) thermal states diagonal in a basis of particle 
or quasiparticle Slater determinants, associated with dif¬ 
ferent values of the running effective order parameters. 
Hence, E'^®p(pspa) = 0, in agreement with the fact that 
PSPA contains just static fluctuations around mean field. 
Such correlated but still unentangled approximation can 
be derived from the auxiliary field path integral repre¬ 
sentation [U, and becomes exact at sufficiently high T 
[13. Its breakdown at low T reflects the onset of entan¬ 
glement, i.e., of a finite value of E‘^®P(p). Eq. ( 1571 ) defines 
a limit temperature above which = Q. Mixtures 
of fermionc gaussian states are also important in noisy 
fermionic quantum computation models (lH . [35l| . 


III. THE CASE OE FOUR SINGLE PARTICLE 
LEVELS 


We will now examine in detail the special case of 
a fermion system with single particle space dimension 
n = 4. This is the lowest dimension where non-trivial 
fermionic entanglement arises, i.e., where E^^p can be 
non-zero, as will be verified. We will extend the results 
of [3, which considered just pure or mixed states with 
a definite fermion number, to general states which do 
not necessarily have a definite fermion numb er, y et still 
have a definite number parity P (see also [^llj]). This 
sp space can accommodate 8 linearly independent pure 
states of the same number parity, so that the Hilbert 
space dimension for fixed P is 8. 
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FIG. 1. (Color online) Schematic representation of pure 
fermion states with odd (top) or even (bottom) number parity. 
A general state with definite number parity is a superposition 
of the 8 states indicated on the left of the dashed vertical line, 
where a full disk indicates an occupied level. In the normal 
form (p'^^P diagonal), obtained after a suitable Bogoliubov 
transformation, it can be reduced to the superposition of two 
states like those indicated on the right. The state is entangled 
(in the sense of not being a quasiparticle vacuum or Slater De¬ 
terminant) iff the product C (Eq. (1431) 1 of the coefficients of 
the left and right groups of four states is non-zero, implying 
nonzero weight for both states of the normal representation. 


A. Pure states 

1. Odd parity states 

We first consider pure states |4f) of this system with 
odd number parity: P|4f) = — Ifk). These states are then 
linear combinations of single fermion states and three- 
fermion states, so a general odd state can be written as 
(Fig. m top) 
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|4/)=^(a4|0)+Ac.|0)), (38) 

i=l 

where |0) = cJc^CgcJlO) is the completely occupied state 
and |Q:p-|-|/3p = 1, with a, j3 four dimensional complex 
vectors. It is easily seen that the single hole states Ci|0) 
are 


^ ^ ^ GjfczcJ-44|0), (39) 

j, k,l 

where Cijki denotes the completely antisymmetric Levi- 
Civita tensor in dimension 4. The elements of the gener¬ 
alized one-body density matrix (EU are then given by 

^ ( 40 ) 

^‘^3 — ^ ^ 5 ( 41 ) 

k, l 

i.e., = aaf — -I- |/3p/4. We now show that the 

ensuing eigenvalues fy of the 8x8 matrix are four- 
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fold degenerate and given by 

,42, 

where C'(|5')) is fully determined by the S 2 entropy (IMl) . 

Cd^-)) = ^S2ip<i^P)/4. = ^tr [p«P(/4-p‘*P)-Kt«;] 

4 

= 2\f3^a\=2\Y,M, (43) 

i=l 

and plays the role of a pure state fermionic concurrence. 
It satisfies 0 < C < 1, and as will be seen in the next 
subsection, it is the generalization of the Slater correla¬ 
tion measure defined in i[ni for two fermion states. 
It also coincides with the quadratic invariant derived in 
using a spinors classification based approach. The 
entanglement entropy (|25]) becomes 

= 4h{f+) = -4(/+ log 2 /+ + /_ log 2 /-). (44) 


Proof. We first consider a unitary transformation c —>■ 
C7c of the operators cj, such that 

a^U^a, /3^Det[C/dC/^/3, (45) 


in (1381) . which does not affect the value of Cd^/)) (Eq. 
(Pl)l. By choosing an orthonormal basis of such 
that the original vectors a and j3 are generated by the 
first two elements (for instance, ei oc a and 62 oc 
/3 — {oL^(3 )ol/\ol\’^), we can use this first transformation 
to set as = 04 = 0, = /34 = 0 in the new basis. In this 

case, Eqs. (|in i) -(FT |) lead to 




/|ad^ + |/32p aid2-/3i/32 
0201-/32/31 |o2p + |/3iP 
0 0 
\ 0 0 


0 0 \ 
0 0 
W 0 

0 I/3IV 


,(46) 


/O 0 0 

0 0 0 

0 0 _ 0 

Vo 0 di/32 — d2/3i 


02/3i 



(47) 


It is then seen that the diagonalization of p'^^P is achieved 
through i) a unitary transformation of the operators ci, 
C 2 , 


Cl = uai + va 2 , C 2 = —vai + ua 2 , (48) 

with |"| = ^ 2 ~(//-^T ^ ^ which 

diagonalizes the first 2 x 2 block of p®P and 1 — /?p, plus 
ii) a Bogoliubov transformation of the operators C 3 , C 4 , 

C 3 = u'oa -I- v'a\ , c\ = —v'as + u'a\ , (49) 


determinant C'^(|4'))/4, leading then to the same eigen¬ 
values f± of Eq. (HU) ( a 2 X 2 matrix with trace t and 

determinant d has eigenvalues *^^*2 □ 

Note from STD that if p‘i®P is diagonal (p®P diagonal 
and At = 0 ) and C'd'k)) < 1 , then necessarily 0,2 = /32 = 0 
or Qfi = /3i = 0 in (SZl). This implies that after the pre¬ 
vious transformations, |'I') can be written in the normal 
form (top right scheme in Eig. [IJ 

14-) =a'al|0a)+/9'ai|5a), (50) 


i.e., P' oc a', with | 0 a) the vacuum of the a operators, 
|0a) = a|444|0a), and \a'\l = /+, = /_ if \a'\ > 

1/3'I, such that C(|'I')) = 2\a'j3'\. This state leads to 



/l«T 

0 

0 



0 

0 

I/3'|"/3 0 

0 |/3'|2 




0 

0 



On the other hand, in the maximally entangled case 
Cd^A)) = 1, /± = 1/2 and p'^^P = Ig,/2 in any basis, 
i.e., after any Bogoliubov transformation. In this case 
P = with |a| = \P\ = l/v^, and the form (l50l) is 

obtained just by choosing ei in the direction of a. 

It is apparent that if /3 = 0 in (l38l) . IdA) can be written 
a single fermion state a|| 0 ), where a| = 
ilarly, if a = 0 I^a) can be written a single hole state 
ai|oV with ai = Accordingly, C(|4')) = 0 in 

these cases. The vanishing of C(|4')) for nonzero but 
orthogonal a and P (Eq. (|43ll l generalizes the previous 
result, showing that in this case |'I') can still be written 
as single quasiparticle (/3' = 0) or quasihole (a' = 0) af¬ 
ter a suitable Bogoliubov transformation of the original 
operators. This includes the three level case where, for 
instance, the fourth level is empty, which implies 04 = 0 
and /3i = 0 for * = 1, 2, 3, leading necessarily to /3 ^q: = 0. 

We also mention that the four eigenvalues of p®P in 
Eq. (|46ll are f± and |/3p, the latter two-fold degenerate. 
Since C'(|4')) < 2|a||/3|, 


/+> 


1 +Vl-4|a|2|/3|2 
2 


Maxjlap, |/3|d , 


being then verified that the eigenvalues of p®P are ma¬ 
jorized by those of p‘^®P and hence, that > S"^®p, 
> S-fP = dC^I^-)). 

Dualization. Eqs. (l38ll and (l43ll indicate that the state 
cdO) plays the role of partner or dual of the state cj|0). 
We may obtain the partner state with the hermitian op¬ 
erator 


T = - 


1 


eijM[c\c]clci + c\cjCkCi] 


(51) 


with |",| = e' = |/3^| — which diago¬ 

nalizes the rest of p^^P, comprising again two 2x2 blocks 
([fl These four 2x2 blocks have all trace 1 and 

'■±K34 \a\^f 


such that for / = 1,..., 4, Tcj|0) = Ci|0), Tci\Q) = cj|0). 
We can then express Eq. 63]) as 

C(|4/)) = |(§|4/)|, |§)=T|vI-)* 


(52) 
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where = '^jOiic\\Q) + ^iCi| 0 ) denotes the “conju¬ 
gated” state in this basis. Note that the 8 x 8 matrix that 
represents T in the basis (c||0),..., C 4 IO), ci|0),..., C 4 IO)) 
is just 


T = 




(53) 


A generalization of (151|) for higher dimensions is consid¬ 
ered in [ 2 ^. 


2. Even parity states 

We now consider pure states of even number parity, 
P|xp) = |xp). They can be obtained, for instance, by 
changing a particle for a hole in the odd-parity states. 
An even state is then a linear combination of the eight 
states shown in the bottom plots of Fig. [I] comprising the 
vacuum | 0 ), six two-fermion states and the completely 
full state |0) = cIcJcgcjjlO). We can write this state as 

4 

!«') = ai|0) -/3i|0) -t-^ajcJ-cllO) -h^jCiCjlO), (54) 

t=2 

which is just Eq. (l38l) with the replacements c\ -O- ci and 
|0) O c||0), implying |0) o — ci|0). Notice that 

cicj |0) = ^ XI 10) ■ (55) 

k,l 

In this notation, the eigenvalues of are then given by 
Eq. (l42l) with the same expression (l4^ for C(|4')). The 
entanglement entropy is given again by Eq. (PI). 
Notice, however, the minus sign in the term associated 
with Pi. The expression (H51) reduces to that of 0 for 
the case of two-fermion states (oi = Pi = 0 ). 

The state ((Ml) is then a Slater Determinant or quasi¬ 
particle vacuum iff Cd^t)) =0. As a check, the quasi¬ 
particle vacuum (EH) corresponds in the present case to 

a oc (1, T 21 , T 31 , T 41 ), 

P oc ( —T 2 ir 43 — T 31 T 24 , — 74ir32, T 43 , T 24 , 732 ) ,(56) 

being verified that — 0 - is also seen that in 

the three-level case (i.e., level 4 empty, implying 04 = 0 
and Pj = 0 for j = 1,2,3) C'(|4')) is always zero. 

The normal form (1501) becomes here 

I4-) =a'|0„)-/3'|0a), (57) 

i.e., a superposition of the vaeuum and the maximally 
oeeupied state (bottom right scheme in Fig. [IJ for the 
diagonalizing quasiparticle operators. Of course, after a 
trivial particle hole transformation Oj -O- a| for j = 1 , 2 , 
we may always rewrite (EH as a sum of two two-fermion 
states, i.e., 


(58) 


which extends the results of @ valid for two-fermion 
states to arbitrary definite parity states. 

The dualization operator m takes here the form 

T = -c\c\c\c\ - C 4 C 2 C 3 C 1 - ^ X ^i 3 kic\c]ckci , (59) 
which satisfies 

T|0) = -|0), r|0) = -|0), Tclc]\0) = ix«bfe;44l0)> 

k,l 

i.e., Tcjc||0) = cicdO), TcicdO) = c|c||0). 

It is represented in the special basis {|0), c^cj |0), 

4 ^ 110 ), 4 c| 10 ), - 10 ), cicj 10 ), 410), 10 )} by the same 

matrix ((501) . We can then write again Cd^t)) in the form 
(150() . If oi = /?! = 0, the ensuing expression reduces to 
that of 0 . 

The two-fermion states considered in i,[ni are only 
a particular case of the more general even states (1501) . 
For two fermion states the contractions (ctCj) obviously 
vanish (k = 0 ), and the eigenvalues fi, of the generalized 
one body density matrix reduce to those of the one 
body density matrix p®P, implying S'^p = S'‘i®p. 


B. Mixed states and analytic evaluation of the 
concurrence 

The fermionic concurrence for mixed states can be de¬ 
fined by the convex roof extension of Eq. ((^ 51 ) . For two- 
fermion states an explicit expression was derived in Q. 
We will here extend this expression to the present general 
states (see also [l^). Let 

p = XAfc|^'fc)(^fe| (60) 

k 

be a mixed state with eigenvectors and eigenvalues 
Afe, with Afc > 0 for fc = 1 ..., r and r < 8 the rank of 
p. We will assume that all |4'fc) have the same number 
parity, such that they are of the form ( 1551 ) or ( 1551 ) . i.e., 
I^fe) = Z)i=i“fci4l0) + PkiCi\0) in the odd case. Ev¬ 
ery convex decomposition p = can be 

obtained from these eigenvectors through a r' x r ma¬ 
trix U with orthonormal columns ([/If/ = A) such that 
y/P]\^]) = UjkV^\^k)- Note that the states |$j) 

are normalized, so that pj = X)/c=i ^kWjk]'^■ 

The average fermionic concurrence (generalized Slater 
measure) of such decomposition is 

{C{{p„ |4>4})) = 

3 3 

= XIX^2-^: ;7\/AfeAi"(^fe|'I'/)| (61) 
3 k,i 

The matrix C of elements 

Cki = 4AfeAi(l'fc|4'i) = 4 XkXppjoik + Piop ), (62) 


|4') = a'44|0a) -I- P'a\a\\Qa) , 





is complex symmetric. Therefore, it admits a decompo¬ 
sition of the form 0 C = VDV"’", where is a unitary 
matrix and Z) is a real diagonal matrix whose diagonal 
elements d* > 0 are the square root of the eigenvalues 
of CCt = CC, sorted in descending order. Defining 
S = UV, Eq. (l 6 T|l then reads 

(63) 

j k 

Since \J2k SjkM > - J2k>2\^jk\dk) = 

di — ^ necessary condition for the “separabil¬ 

ity” of p, i.e., for p to be a convex mixture of Slater 
determinants with the same number parity, is 

di <'^dk. (64) 

k>2 

As in the case of two fermion states, we will now show, 
following the scheme of Q, that this is also a sufficient 
condition for separability. Indeed, from (1551) it is seen 
that p is separable if there is a matrix S with orthonormal 
columns such that for every j, 

r 

\Y,dkS]k\=0- (65) 

k=l 

Now, provided condition (l64l) is fulfilled, there are always 
phases 9k, k = 2 ,..,r such that di = | X)fc= 2 I- 

Then a matrix with elements Sjk = -—i where 
Pjk = 0,1 and 9i — 0 , will give the desired decompo¬ 
sition if the signs can be arranged such that the 

condition S^S = Ir is satisfied. This can be ensured 
by taking r' = 2 if r = 2, r' = 4 if r = 3,4 Q and 
r' = 8 if 5 < r < 8 , where we can set pji = 0 V j and 
(/iifc,...,/i 8 fe) as ( 0 , 0 , 0 , 0 , 1 , 1 , 1 , 1 ),( 0 , 0 , 1 , 1 , 0 , 0 , 1 , 1 ), 
( 0 , 0 , 1 , 1 , 1 , 1 , 0 , 0 ),( 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 ),( 0 , 1 , 0 , 1 , 1 , 0 , 1 , 0 ), 
( 0 , 1 , 1 , 0 , 0 , 1 , 1 , 0 ),( 0 , 1 , 1 , 0 , 1 , 0 , 0 , 1 ) for k = 2 ,..., 8 . 
This completes the proof. 

On the other hand, if condition (IMl) does not hold, the 
average (1631) is not smaller than di — X)fc =2 ■ This lower 
bound may be achieved with the same construction used 
above, choosing 6 k = 7 r /2 for k > 2. Then, the mini¬ 
mizing decomposition is that where all the components 
have the same concurrence, which is the concurrence of 
the state p, 

r 

C{p) = Min{p^.j$.)}EPiC'(l^j)) = Ma.x[di-'^dk,0] ■ 

j k=2 

( 66 ) 

Using the dualization matrix (1531) we may also obtain the 
eigenvalues dk as those of 

R = v'pi/2Tp*Tpi/2 , (67) 

where p* means conjugation in the basis where T takes 
the form (1551) . 


Once C is obtained, we can evaluate the convex roof 
extension (1571) of as 

£;'i"P(p) = ’ (68) 

in the same way as in the two-qubit case (^ . since for 

pure states we have similarly = 4/i.( ^~*~^^ ^ (|i')) ^ 
(Eq. (1441) 1. which is a convex increasing function of 

C'(l'k)). The quantity ~ is also the maxi¬ 

mum fidelity between p and a convex mixture of gaussian 
states, as shown in | 22 l| with a different treatment based 
on group-theoretical methods. 

A general mixed state p satisfying [p, P] = 0 will be a 
convex mixture of pure states with even and odd number 
parity. It can be written as a convex mixture of even and 
odd parts, i.e., 

P = V+P++V-P- , (69) 

where p± = 2 ^(i ^ even and odd com¬ 

ponents of p and p± = Tr p(l ± P)/2 the corresponding 
probabilities. Since we just consider pure states with def¬ 
inite number parity, for the general mixed states (j69l) we 
may just take P'i"P(p) = p+P^"P(p+)-kp-P"'"P(p-), with 
pa^P^pj.) evaluated with Eqs. (l6^ and (l68l) . 

As illustration, we consider a definite parity mixture of 
a maximally entangled state jdt) (C(|'I')) = 1) with the 
fully mixed state, 

p=p\^){^ + {l-p)h/B, (70) 

where 0 < p < 1. In the odd parity case, |'I') can be 
written in the form 

14/) = ^(cl|0) + ci|0)) = ^(cl + clclclm , (71) 

whereas in the even parity case we can take |4/) = 
-^(|0) -I- |0)) or -(^( 4^2 + 44)16)- A direct calculation 
using (ESI) leads to 

C(p)=Max[V.0], (72) 

indicating entanglement for p > 3/7, i.e. q > 1/2, 
where q = (4/|p|4/) = p -|- (1 — p)/8 is the total weight 
of 14/). A similar calculation but considering just two- 
fermion states, p 2 = p|4/)(4/| -I- (1 — p)Iq/6, leads in¬ 
stead to C(p 2 ) = Max[ E implying entanglement 
above a slightly smaller value of p (p > 2/5, entailing 
again q = p + [l — p)/Q > 1/2), with C{p 2 ) > C{p) for 
p G (2/5,1). As in the two-qubit case, the existence of 
a finite threshold probability p for non-zero C and hence 
pqsp, implies a finite limit temperature for entanglement 
Tl if p represents a thermal state ( (i-p)/8 ^ 
with Eq the energy of |4/) and Ei > Eq that of remaining 
7 levels), which is larger in the second canonical case. 
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IV. CONCLUSIONS 


We have presented a general consistent formalism 
for describing entanglement-like correlations in general 
fermion states with no definite fermion number yet fixed 
number parity. We have first defined a single level en¬ 
tanglement entropy that quantifies the entanglement be¬ 
tween a single-particle mode and its orthogonal comple¬ 
ment, through the definition of suitable reduced states 
for such a partition of a given basis of the single-particle 
space. The sum over all sp modes of this entropy, Sc, 
can be taken as a measure of the total entanglement of 
the system with respect to this basis, and its minimum 
over all sp bases, was shown to be a function of 
the one-body density matrix, being then invariant with 
respect to unitary transformations in the single-particle 
space. Moreover, if minimization is extended over all 
quasiparticle basis, the resulting entanglement entropy, 
is a function of the generalized one-body density 
matrix, remaining therefore invariant under general Bo- 
goliubov transformations. Such entropy vanishes iff there 
is a single particle or quasiparticle basis in which ev¬ 
ery level is separable from its orthogonal complement, 
i.e., iff each of these levels is either empty or occupied. 
These entanglement entropies satisfy the inequality chain 
Sc > S^P > S^^P. The convex roof extension of S^^p was 
also introduced, its vanishing rigorously identifying “clas¬ 
sically” correlated mixed fermion states which can be ex¬ 
pressed as convex mixtures of pure states or quasiparticle 
vacua, like those emerging at sufficiently high tempera¬ 
tures in interacting many-fermion systems through ap¬ 
proaches like the SPA. 

In the case of fermion systems with four single particle 
levels, a fermionic analog of the two-qubit pure state con¬ 
currence was defined in terms of p^^P, which reduces to 
the Slater correlation measure defined in [iini for two- 
fermion states. The eigenvalues of the generalized one- 
body density matrix, which are four-fold degenerate, can 
be written as functions of this concurrence and conse¬ 
quently, the entanglement entropy S'^^^p is related to the 


fermionic concurrence by an expression analogous to that 
of the two-qubit case. This result suggests that “particle 
entanglement” may be seen as a minimum “mode en¬ 
tanglement” . For mixed states with fixed number parity 
of this system, an explicit expression for the fermionic 
concurrence, defined as the convex roof extension of the 
pure state concurrence, was derived, in complete anal¬ 
ogy to the two-qubit case, which generalizes the result of 
and provides a closed analytic expression for the 
convex roof extension of S'‘^®p. 
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Appendix A: Quasiparticle vacuum 

According to Thouless theorem the vacuum |0o) 
of the quasiparticle fermion operators m is given, if 
Det t/ 7 ^ 0, by [H 

|0a) = 7exp[i T,jc\c]]\0) 

= 7 [1 + 5 T^jclc] -h .. .]|0), (Al) 

where 7 = y|DetC/| and T = —U~^V is an antisymmet¬ 
ric matrix, with |0) the vacuum of the Cj operators. Eq. 
(El) can be verified by directly applying to (|Aip (if 
Det 17 = 0, |0a) can be obtained by applying additional 
creation operators cj to Eq. m)- 

If l^f) = |0a), then /i, = (0a|aj,a,.|0a) = 0 V imply¬ 
ing S’^^P = 0. However, it is easy to see that 

p^P = I-(0jcct|0„) =UI/t^ (A2) 

implying S^p > 0 if U 7 ^ 0. The eigenvalues pk of p^P 
are then just the square of the singular values of V. The 
state |d>) appears, therefore, mixed at the sp level, re¬ 
flecting that it cannot be written as a Slater determinant 
in operators of the form (IT^ . 
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